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Abstract 

The linear response conductance coefficients are calculated in the scattering 
approach at finite frequency, damping and magnetic field for a microstructure 
in which the reservoirs are modeled as quantum wire leads of infinite length 
but finite width. Independently of frequency, inelastic scattering causes sub- 
bands with large group velocity to contribute more strongly to the conduc- 
tance than channels of comparable transmission but slower propagation. At 
finite frequency and magnetic fields, additional correction terms appear, some 
of which are sensitive to the phase of the S matrix. 
The Landauer-Biittiker formula [0] for the low-temperature, linear-response DC conduc- 
tance coefficients of a small conductor requires only the knowledge of the quantum mechan- 
ical transmission coefficients, i.e. absolute values of S-matrix elements. The S-matrix in 
turn appears because one necessarily has to consider the system as being open, i.e. having 
a continuous spectrum. With discrete levels, no dissipative response is possible. One often 
uses a model in which the sample is connected to the outside world by straight quantum 
wires of infinte length 0]. 

Here, we want to study the linear conductance of such an open system under an AC 
perturbation. The limiting case of infinitely broad reservoirs and no damping has been 
considered previously, , yielding an admittance matrix ioi n ^ I 



Here, n and / label the leads so that = J2 Qniii^) Vi(ct;), and the trace extends over 

I 

subbands. An apparent conflit has been noted |^ between the widespread use of Landauer- 
type formulas even in systems that are one- dimensional at infinity, and the requirement that 
the reservoirs have to be much wider than the sample itself to permit equilibration. 

We provide an answer to this open question by considering the linear response of a system 
with quantum wire leads, including an inelastic scattering rate in the equation of motion. 
Electron-electron interactions are assumed to be absent in the reservoir region. An external 
uniform magnetic field B is taken into account. 

Consider a system under the infiuence of an external perturbation, H = Hq + Hi, where 
Hq is the unperturbed time-independent Hamiltonian and (in the interaction picture) 

Hu{t) = ~ej n/(x, t) <l>(x) rfx = e'^' Fi{t). (2) 
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Here, nj{t) is the number density operator. The assumption that a static potential $ can 
describe the field distribution amounts to neglecting the rotational part of the time-varying 
electric field. The fact that the resulting perturbation (which must also contain the self- 
consistent induced fields) just adds to the original Hamiltonian Hq is a consequence of the 
linear response approximation. The expectation value of the current density is (j(x, t)) = 
Tr pit) j, where we use the linear approximation in the von-Neumann equation for the density 
matrix, 

^^pi{t) = 1 [Hu{t),po]-l {pi{t)-po}. (3) 

Here Pq is the equilibrium density and 7 is an inelastic relaxation rate . To solve this in 
the basis of many-body eigenstates |a) of Hq with energy e^, one needs the matrix elements 
of Eq. (H). Using the continuity equation, hi{-x,t) = —V ■ j/(x, t), one can perform an 
integration by parts as in Ref. [0], yielding 

{a\Fj{tm = I / (a|j,(x,t)|/3)$(x)rfx 

- j <l>(x)V-(a|j/(x,t)|/3)rfx|. (4) 

Here, j/ is the current operator associated with Hq, containing at most a static vector 
potential if the Landau gauge is chosen in each lead 

With this procedure, we have traded the infinite integration domain of Eq. (|]) for a 
finite domain A. The price that was paid is the additional boundary integral over dA. At 
u; = 7 = 0, the contribution of the bulk term is known to vanish p|. Its significance at 
finite frequency as compared to the surface term is not addressed in Ref. , but it has been 
stated H] that it should be negligible at small frequencies, when one is in the quasi-stationary 
regime, Vj ~ 0. However, an estimate of the conditions for this regime is not available in 
the literature. 

Using (q;|p|/3), the current response becomes 

t)) = jo(x) - e 2^ — - — ■ — [Pa-Pp). (5) 

^ hu + ea - + ih'-f 

Here, jo(x) is the unperturbed equilibrium current, and Pa = (a|po|«)- The total current 
is the integral of Eq. (|^) over the cross section of lead n. Furthermore, we only need the 
deviation from the equilibrium value, so that we bring the term involving jo(x) to the 
lefthand side. The resulting induced current has the time dependence e*"^*, which we drop 
to get the corresponding Fourier component (which can still be a function of the distance x 
from the sample). 

For the bulk term in Eq. (B), we can calculate the total current through all the leads. 



E/r' = / v-(r^x')-jr(x'))rfx' 

A 
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J {(3\n{x.')\a) $(x) (a|n(x)|/3) dx'dx, (6) 

A 

where Gauss' theorem and the continuity equation were used. We henceforth abbreviate 
z = uj — i'-f. To hnear order in z, this can be expressed in terms of the Green's function, 
= E«|a)(«|/(^-ea), as 

^Jbuik = _,e2^ Y.Pfs [ H^) (/3|r2(x')G'(e^)n(x) + h.c. |/3)dx'dx. (7) 

The similarity of this operator expression to the mean-field-hke result of Ref. [cf. Eq. (54) 
therein, summed over leads] leads us to identify this contribution with the internal response, 
which we have thus re-derived in an alternative way. 

In the following we shall focus on the boundary term and find a generalization of the 
external response result of Ref. [^0. The mutually decoupled, straight infinite quantum 
wires leads at large distance from the sample are described by a non-interacting Hamiltonian 
Hqw The many-body eigenstates |a) of the full unperturbed system decompose into Slater 
determinants of single-particle eigenstates of Hqw in the asymptotic region, a complete set 
of which is given by the scattering states, ip^ap at energy E, consisting of an incoming wave 
in subband a of lead p, and of outgoing waves in all the subbands of all the leads /. Using 
local coordinates x, y parallel and transverse to the wire /, the energy-normalized quantum 
wire eigenfunctions are 



da(x) 



2tt 



dE: 



dk 



-1/2 



xkJy)- (8) 



Unless i? = in the leads, the transverse wavefunctions xiaiv) are also functions of the 
wavenumber k^^(£'). The + (— ) sign denotes propagation toward (away from) the sample. 
The scattering states in lead / then are 

^Eap = 5pl da(x) + E Slp,a'a ^U' • (9) 
a' 

Here, Sip^a'a is the 5* matrix element for scattering from subband a of lead p to subband a' 
of lead /. Bound states of the system do not enter in the completeness relation because of 
their exponential decay. 

After making the transition to the single-particle scattering states, the boundary contri- 
bution to the current at a distance x along lead n from the sample takes the form 



N N 



Ux) = ~te'hY,ViY. Yl jdE'dE 



1=1 P,q=l a,b=0 , 



f{E') - f{E) 
E-E' 



E'-E-hz • ^^"^ 
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Here, is the component of the single-particle current operator along lead /, integrated 
over the cross section at a distance Xyi from the sample (determined by OA). We use the 
convention that subscripts a or a' are always associated with a dependence on while b 
or h' labels functions of E' . If we now insert Eq. (^), then the last line in Eq. (|10D contains 
a product of current matrix elements between quantum wire eigenfunctions at different 
energies, of the type considered in Ref. 0. For brevity, we restrict ourselves to Qni with 
n ^ I. In that case one can exploit the 5pi in Eq. (^, as well as the unitarity relation 
S'^{E) S{E) = 1. Furthermore, current matrix elements between counter-propagating waves 
such as {^'_a\J''\^+b) vanish as a; — > 0. To linear order in a;, this leaves 



gni = -ie^h / dE' dE 

1. Ill r\ J 



fiE') - f{E) 



a,b,a',b'=0 , 



E-E' 

E' — E — hz {^^-"I'^'l^-'''^ {^+b'\J"'\^+a') Snlfi'bSnl,a'a 
+ {^+a'\j''\^+b') {^-blJ^l^-a) Sln,a'aSln,b'b] (H) 

The arguments on the current operators were left out, keeping in mind that the coordinates 
are in leads / and x in lead n. Following Ref. 0, one of the energy integrals is done such 
that all terms multiplied by f{E) are integrated over E', vice versa for terms multiplied by 
f{E'). Applying the residue theorem, the poles of the Fermi function then do not enter. If 
the first term in braces is integrated over E {E'), we close the contour by a large semicircle in 
the positive (negative) imaginary plane; for the second term in braces, the opposite contour 
is taken. With this choice, all the e^^^^ factors in the wire eigenfunctions ^ yield exponential 
suppression along the return contour due to the imaginary part of k. Poles of S occur 
only in the negative energy plane, and in the positive plane for S*. As a consequence, the 
above choice of contours always selects the half plane in which no poles of the S-matrix are 
enclosed. The contributions from the pole at E' = E cancel due to the difference of Fermi 
functions, so that the only pole that remains is the one ai E' = E + hz. However, this pole 
is never enclosed for the second term in braces, so that the latter makes no contribution at 
all. We are thus left with 



oo 

gni = 2ne^h ^ dE 

1.11/ n 



a,b,a'b'=0. 



f{E + hz)-f{E) 

hz 



'X {^~a\-J^\^~b) {^lb'\-J''\^+a') S*i^f,,f,Snl,a'a ■ (12) 

At i? = 0, the transverse modes are orthogonal irrespective of energy, yielding immediately 
the requirement a = b, a' = b'. Consider these terms first, but at i? 7^ 0. Expand all 
wavenumbers around E, recalling that b labels functions of E + hz, e.g. 

kUE + hz) ^ kUE) + qi,, gi^^^^. (13) 



v: 



-aiEY 



Here, we introduced the channel- specific group velocity Vg _^^{E) = {1/h) dEl^{k}^^{E)) / dk. 
Then first-order perturbation theory yields for the transverse channel eigenfunctions 

X-b,E^n. - X-a nu^cq-a ^^^X-a'^. _ El.ikUE))' ^'^^ 
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where ujc = eB/mc. This can be used to evaluate the current matrix elements in Eq. 
(|1|) with the Landau gauge, where one needs {xia\X-b), {X-alv'lX-b) and (x+6'lx+a')> 
{x+b'\y\x+a') ■ To get explicit expressions, we specialize to parabolic quantum wires with 
the same dispersion relation in leads n and /, 



Ea{k) 



27.2 



2m cUq + uj^ 

The terms in Eq. (|12D with a = b, a' = b' then yield to linear order in z 

f{E + hz)-f{Ey 



diag 
9nl 



e 
J 

( 

xe 



E 



dE 



hz 



SnLa'aiE + hz) Snl,a'aiE) 



(15) 



(16) 



[x/\v^J+x^/\v'gJ 



TTlLUt. 
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This reduces to Eq. ([l|) only if both Uc = and one drops the exponential, which describes 
damped oscillations as a function of x and x_4 since z = — i 7. Note that Xj^ 7 is the 
velocity that a particle must have in order to traverse the distance x_4 before the inelastic 
damping becomes significant. Similarly, Xj^u is the velocity that is required to traverse x_4 
in one oscillation period of the external field. If carriers can then enter and leave the sample 
region so fast that zx/\v\ 



ga' 



zxj\^/\Vg^ ~ 0, then the precise location of the current and 



voltage probes, which is given by x and x_4, respectively, becomes irrelevant. In particular, 
channels exceeding the cutoff group velocity, max(a;7, Xj^'j), are essentially undamped in Eq. 

Their number increases when the reservoirs are much wider than the sample region, in 
which case "slow" subbands typically have negligible transmission through the sample since 
coupling to the narrow sample requires large transverse momentum transfer. For narrow 
reservoirs, the transmission of such slow subbands can be appreciable, so that the deviation 
between Eqs. (|I]) and (|T^ can be significant even at a; = i? = 0. 

We also find contributions from off-diagonal terms with a ^ b or a' ^ b' in Eq. (|12|), of 
which we list only one example [ITO]: 



off 
9nl 



h uJq 



2 

e uJc 



I h m ^ 
2mujo h 



dE 



f{E + hz)-f{E) 



hz 



X 



nl\a' ,a+l 



{E + hz) SniME) 



(17) 



This contains further oscillatory exponentials, which for large Xj^ lead to suppression of the 
integral. However, this may be compensated by the fact that g!^{^ has roughly one power 
of wavenumber more than g'^i"'^. In that case the admittance at B ^ can give information 
about products of S-matrix elements that cannot be written as a trace over subbands as in 
Eq. (0). 

We benefited from discussions with M. Biittiker, M. Leadbeater, E. McCann and 
A. D. Stone. 
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